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ABSTRACT 

In real telescopes, the optical parameters evolve with time, and the degradation is 
often not uniform. This introduces variations in the image profile and therefore 
photo-centre displacements which, unless corrected, may result in astrometric 
errors. The effects induced on individual telescopes and interferometric arrays 
are derived by numerical implementation of a range of cases. The results are 
evaluated with respect to the potential impact on the most relevant experiments 
for high precision astrometry in the near future, i.e. Gaia, PRIMA and SIM, and 
to mitigation techniques applicable from design stage to calibrations. 
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1 INTRODUCTION 

The mirrors of ground based instruments suffer degrada- 
tion of the reflecting surface, e.g. by oxidation or contam- 
ination of the coating, due to dust, moisture and chemical 
agents in the environments (jVucina et al.|[2006l ). The re- 
sult is a degradation of the optical throughput, which in 
general is not uniform, but rather described by a com- 
plex, often patchy, distribution. Also, the overall tele- 
scope transmission is progressively reduced, so that it 
becomes necessary to process in particular the primary 
mirror regularly for removal of the degraded reflecting 
coating and application of a new one. 

Space optics is expected to be much less affected by 
surface characteristics variations, thanks to the stable en- 
vironment and suitable protective layers, but degradation 
e.g. due to 7 radiation is still experienced, up to 1% of the 
overall throughput, d epending on the choi ce of substrate 
and coating material l|Baccaro et al.|[2004 ). 

Even small variations are still potentially relevant 
with respect to high precision astrometric measurements, 
and as such their effects have been investigated and are 
discussed in this document. The main experiments in 
which astrometric errors induced by non-uniform optical 
throughput variation are potentially relevant are Gaia, 
PRIMA a nd SIM. 

Gaia jPerrvman et al.ll200ll ) is the most important 
European space mission devoted to micro-arcsec {fias) 
astrometry currently being implemented. Interferometry 
also aims at comparable levels of precision by means 
of both ground based arrays, e.g. the Phase Referenced 
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Ima ging and Microarcsecon d Astrometry [PRIMA] facil- 
ity (|Delplancke et al]l2003l ) of the Very Large Telescope 
Interferometer [VLTI] of the European Southern Obser- 
vatory [ESO], o r from space, e.g. th e Space Interferometer 
Mission [SIM] (|Unwin et al.ll2008l ). 

The astrometric performance of an astronomical in- 
strument depends on the detected signal profile. The 
Point Spread Function (PSF), i.e. the intensity dis- 
tribution from a point-like object at infinity, is built 
from the d iffraction integral, de scribed in optics text- 
books, e.g. iBorn and Wolj (|l985l ): the numerical imple- 
mentation used here was described in a previous paper 
(IGai and Cancelliere ''2007^. The PSF includes, by a suit- 
able wavefront error (WFE) map, the realistic descrip- 
tion of several optical characteristics. The polychromatic 
PSF is produced by superposition, weighted by spectral 
distribution, of the monochromatic PSFs. 

Hereafter, Im is the monochromatic PSF at a given 
wavelength A, and Em is the corresponding complex am- 
plitude distribution on the focal plane; F is the effective 
focal length of the telescope. The coordinates on pupil 
and focal plane are respectively {(,, rj} and {x, y} (lin- 
ear units). The complex amplitude depends on the pupil 
function P and on the aperture transmission function T: 



-iTv(x^ + yr])/XF 



(1) 



The constant k provides the appropriate photometric re- 
sult associated with the source emission, exposure time 
and collecting area. The monochromatic PSF can be ex- 
pressed as 



{x,y) = \Em {x,y)f . 
The pupil function P (^, rj) 



(2) 

depends on the 
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phase aberration function $, often expanded e.g. in terms 
of the Zernike functions ijBorn and WoHlll985l l: 



(3) 



The aperture transmission function T is often considered 
to take unity value inside the pupil and zero outside; it is 
therefore used simply as a convenient numerical tool for 
insertion of the pupil geometric description. In the cur- 
rent simulation, however, it is used also to describe the 
variation of optical throughput of the telescope through- 
ou t the pupil. Th e formalism is similar to that adopted 
bv lLinfieldl (|2004l ) . in the framework of straylight analysis 
for a Terrestrial Planet Finder (TPF) coronagraph. The 
most relevant changes are a simplified treatment, since we 
consider only transmission variations and neglect phase 
contributions, and the focus on astrometric effects, i.e. 
image photo-centre displacements. 

We expect that a uniform variation of the optical 
throughput over the whole telescope aperture does not 
introduce any variation on the image photo-centre; how- 
ever, for a generic aberrated configuration, a non-uniform 
pupil transmission introduces an amplitude modulation 
of the pupil function, thus modifying the resulting fo- 
cal plane energy distribution, hence the measured signal 
profile, and consequently its photo-centre estimated po- 
sition. Such photo-centre displacement is associated with 
an error on the astrometric measurement. Although in 
general the optical throughput is wavelength dependent, 
our analysis is referred to a single spectral type source, 
and uniform spectral dependence is assumed. Further lev- 
els of detail may be introduced in future investigations. 

In section [2] we discuss the introduction of transmis- 
sion non-uniformity features and evaluate a set of ref- 
erence cases focused on the astrometric performance of 
one telescope. In section[3]we analyse some consequences 
on the measurements of an interferometric array. Then, 
in section \^ we discuss the implications of our findings 
on design and calibration of astrometric instruments. Fi- 
nally, we draw our conclusions. 
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Figure 1. Throughput distribution (in percent) for a selected 
patch position. The transmission is reduced locally by ^1%. 
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Figure 2. Histogram of throughput (in percent) vs. patch 
position ranging over the telescope pupil. 



representative of the Gaia nominal configuration and op- 
eration. 

The optical throughput perturbations described be- 
low are applied to a range of cases with different im- 
age quality, achieved by selection of the 21 lowest or- 
der Zernike coefficients, corresponding to the normalised 
radius ra ised up to the fifth power, as in the table in- 
cluded in lGai and Cancellierel (|2005l ). The optical images 
are generated as pseudo-polychromatic, using a central 
wavelength of 617 nm, with weight 100%, and two side 
wavelengths at ±100 nm, with weight 50%, approximat- 
ing the instrument response to a blackbody source at 
r = 3 X 10* K. The photo-centre estimation is performed 
using the centre of gravity algorithm, as it is simple and 
independent of signal models. 



2 SINGLE TELESCOPE 

The analysis is focused on a configuration similar to 
the Gaia telescope optics and operation, also because 
its rectangular geometry is marginally more convenient 
for numerical simulation. The primary mirror size is 
1.45x0.45 m; the effective focal length is EFL = 35m, re- 
sulting in an optical scale of 5.89 arcsec/mm. For simplic- 
ity, the study case in this document is one-dimensional; 
this is also most immediately related to the Gaia elemen- 
tary measurements. The PSF is integrated in the direc- 
tion not used for the astrometric measurement by bin- 
ning on the CCD detector; this is simulated by numeri- 
cal integration over an appropriate region. Also, some of 
the effects associated with the instrument characteristics 
and operations are included in the form of equivalent ef- 
fects on the PSF through the corresponding Modulation 
Transfer Function (MTF); in particular, this includes the 
detector pixel (size 10 x 30 ^m) and Time Delay Integra- 
tion (TDI) readout l|Gai et al.ll 19981 ). The values used are 



2.1 Opacity patches 

The telescope throughput degradation is introduced as 
a localised reduction with respect to the nominal unity 
value (To = 1), i.e. an opacity patch p, so that the lo- 
cal transmission over the pupil is t (^, rj) = To — p (^, rj). 
A Gaussian shape of the opacity patch is adopted, with 
given characteristic size a (a = 5 cm for most of the cases 
below), peak position {^o,'7o} and peak value po: 

P (C, '?) = Po • exp -^^ . (4) 

For all cases investigated, the peak value is set to po = 
1%, i.e. the perturbed telescope transmission ranges be- 
tween 99% and 100% of the unperturbed value. The opac- 
ity distribution is therefore not normalised, and its area 
is Ap — . A selected case corresponding to patch po- 
sition close to the centre of the pupil is shown in figure [T] 
The overall optical throughput variation can be expected 
to be dependent on the peak and size of the opacity dis- 
tribution, i.e. of order of po ■ Ap/Ar, where At = L^^Lr, 
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is the primary mirror area. For a generic centre position, 
part of the patch is outside the pupil, further reducing the 
throughput degradation. The overall transmission varia- 
tion, averaged over the pupil, is thus well below 0.1%, 
because the opacity patch is much smaller than the pri- 
mary mirror size. Thus, photometric measurements are 
affected by marginal degradation. 

2.2 Simulation 

The position of the opacity patch is translated uniformly 
over 100 position on the long side of the telescope, and 
over 10 positions on the short side; each case is thus evalu- 
ated on a total number of 1000 instances. The pupil and 
focal plane sampling resolution is respectively 2x2 cm 
and 2x6 ^m. An histogram of the transmission distri- 
bution is shown in figure (2] The throughput mean value 
over the sample of patch positions is 0.99982, with stan- 
dard deviation 3 x 10~^. Hereafter, a range of WFE cases 
is evaluated against the set of throughput perturbations 
described above. In all cases, the readout area is set to 
12 x 12 pixels. 

2.2.1 Case 1 

A set of random values of the 21 lowest order Zernike co- 
efficients is selected, with normal distribution and peak 
value of 50 nm. The corresponding WFE has RMS value 
of 35.05 nm, corresponding to a fairly good optical qual- 
ity, ~ A/20 at A = 700 nm. The WFE distribution over 
the pupil is shown in figure [3] (top). The photo-centre 
computed on each instance of opacity patch position is af- 
fected by displacements with respect to the unperturbed 
case; its distribution vs. opacity patch position is shown 
in figure [3] (mid). The photo-centre displacement mean 
and RMS are respectively 0.03 and 3.34 ^as; the peak- 
to-valley (PTV) is 14.13 ^as. 

2.2.2 Case 2 

In order to evaluate the dependence on the amount of 
aberration, the WFE shape is retained, but its ampli- 
tude is reduced by a factor two by halving all Zernike 
coefficients. The resulting RMS WFE is thus reduced to 
17.53 nm; the map is not shown here, since with respect 
to figure [3] (top) only the scale changes, by a factor two. 
The photo-centre distribution is similar to that from case 
1 (figure [3] mid), but it is also reduced by a factor two; 
the mean and RMS are respectively 0.02 and 1.67 /las; 
the PTV is 7.09 /^as. The photo-centre discrepancy be- 
tween case 1 and case 2, after multiplication of the latter 
by a factor two, have mean and RMS respectively zero 
and 0.03 ^as. 

2.2.3 Case 3 

Similarly, in case of larger aberration with the same shape 
(achieved by multiplication by a factor two of the Zernike 
coefficients, resulting in RMW WFE = 70.11 nm), the 
photo-centre displacement is doubled, but again simi- 
lar to that of case 1 (figure [3l mid). The mean and 
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Figure 3. Map of WFE for case 1 (top, in nm); distribution of 
photo-centre (in /las) vs. opacity patch position, case 1 (mid) 
and 4 (bottom, aberration reversed in sign). 
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Figure 4. Realistic WFE selected for case 6 (top) and the 
corresponding distribution of photo-centre (in /.tas) vs. patch 
position (bottom). 



RMS photo-centre displacement are respectively 0.07 and 
6.67 /xas; the PTV is 27.92 /ias. The photo-centre discrep- 
ancy between case 1 and case 3, after division of the latter 
by a factor two, have mean and RMS respectively zero 
and 0.09 /xas. 



2.2.4 Case 4 

In case of symmetric aberration, i.e. changing the sign of 
the Zernike coefficients, the same RMS WFE is retained, 
and the focal plane image is symmetric with respect to 
case 1. Correspondingly, the photo-centre displacement, 
shown in figure [3] (bottom), is opposite to the initial dis- 
tribution of case 1: the mean and RMS are respectively 
-0.03 and 3.34 /xas; the PTV is 14.13 ^as. The photo- 
centre discrepancy between case 1 and case 4, after sign 
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Figure 5. Photo-centre distribution (in ^as) of case 6 with 
respect to case 1. Different WFE distributions induce quite 
different astrometric errors. 



reversal of the latter, have both mean and RMS equal to 
zero within 10~^ /xas. 



2.2.5 Case 5 

The diffraction limit case is also worth considering, 
because the throughput non-uniformity introduces an 
asymmetry in the system, which is no longer ideal. This 
potentially might result in a photo-centre displacement. 
The simulation shows that this is not the case, since the 
photo-centre mean, RMS and PTV values are all zero. 
This is consistent with the mathematical framework of 
the diffraction integral, since replacing in Equation [T] the 
null aberration <1> = 0, a unity pupil function is achieved, 
and we can see that sign reversal on the focal plane co- 
ordinates only introduces a global phase inversion on the 
complex argument, which vanishes in the square modulus 
providing an invariant intensity distribution (Equation 
[2}. Besides, in mathematical terms, the Fourier trans- 
form F{lj) of a real function f{t) is hermitian, so that 
F{—uj) = F*{lj) (complex conjugate), with the same re- 
sult. 



2.2.6 Case 6 

In order to test a different shape of aberration, an alter- 
nate set of Zernike coefficients is generated, in the same 
range of values, with a resulting RMS WFE of 35.23 nm. 
The same opacity patch distribution as above is used. 
The WFE distribution is shown in figure 3] (top); it is 
equivalent in terms of RMS value, but different in shape, 
to that of case 1 (figure [3l top). The photo-centre dis- 
placement distribution, shown in figure 3] (bottom), has 
not only different shape with respect to that derived from 
the previous WFE distribution, but also larger deviations 
from the unperturbed case: the mean and RMS are re- 
spectively 0.13 and 5.45 /uas; the PTV is 21.14 /xas. The 
plot of photo-centres of cases 1 and 6 against each other, 
shown in figure [S] in which each point is affected by the 
same opacity patch, does not evidence significant correla- 
tions. The RMS WFE is therefore not in itself a sufficient 
indication of the instrument sensitivity to non-uniform 
throughput variations. 
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Figure 6. Distribution of photo-centre (top, in /^as) vs. patch 
position, case 7 (smaller patch), and photo-centre distribution 
(bottom, in /^as) of case 7 with respect to case 1. 



2.2.7 Case 7 

In this case, the characteristic size of the Gaussian dis- 
tribution of opacity is reduced by a factor two, to cr = 
2.5 cm. Since the peak value remains set to 1%, the effec- 
tive perturbation to the system is smaller, and a reduced 
astrometric effect can be expected. The photo-centre dis- 
placement distribution, shown in figure [S] (top), is con- 
sistent with this expectation: the mean and RMS are re- 
spectively 0.00 and 0.93 ^ias; the PTV is 4.63 pas. The 
photo-centre plot of case 7 vs. 1, shown in figure [H] (bot- 
tom), does evidence a significant correlation, since for 
each point the two cases have in common both the opacity 
patch and the WFE map. Scaling the former by a factor 
3.56, defined by best fit, the discrepancy mean and RMS 
are respectively 0.02 and 0.29 /xas. The scaling factor is 
close, but not equal, to the ratio of the opacity patch ar- 
eas (i.e. four) on the pupil plane; given the non-linearity 
of the diffraction integral, Eqs. [1] and [2j this seems to be 
reasonable. 



2.3 Comparison 

The various WFE cases are associated with different 
photo-centre displacements, induced by the same set of 
telescope throughput perturbations. A compact represen- 
tation of the relationship among the first four cases is 
shown in figure [7] The image profile is more aberrated 
in case 3 than in cases 1 and 2, and correspondingly a 
larger fraction of the total flux falls outside the readout 
region; this introduces additional noise. The astromet- 
ric error induced by a given transmission perturbation, if 
the WFE distribution is uniformly scaled by a factor 0.5, 
2, and -1, is with good approximation multiplied by the 
same amount. 

Independent WFE distributions induce quite differ- 
ent astrometric errors, with respect to the same transmis- 
sion perturbation, as shown in figure[5l Thus, the photo- 
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Figure 7. Photo-centre distribution (in fias) of cases 2, 3 and 
4 with respect to case 1. Scahng the WFE ampHtude induces 
a proportional change in the resulting astrometric error. 



centre displacement depends on the focal plane position. 



Let us consider the case of a two telescope inter- 
ferometer with baseline B, observing a pair of stars with 
separation (p « 1 rad and zenithal distance 9 (for the ob- 
ject with higher elevation) , in a planar geometry for con- 
venience. The variation of optical path difference (OPD) 
between the two stars, measured by the interferometer, 
is then AOPD ~ ip B cos9, so that the baseline error 
i7{B) — (t(^) is reflected into an astrometric error cT{ip) 
on the determination of the star separation ip: 



B 



(8) 



Such effects are relevant e.g. for orbit determination, 
when measurements taken at different epochs ti, t2 are 
combined; in the mean time, the instrument may have 
suffered changes not easily identified (e.g. < 1% in op- 
tical throughput) but sufficient to introduce relevant as- 
trometric errors. We can evaluate the magnitude of the 
effect on two cases, namely PRfMA at the ESO VLTI 
(ref. PRIMA) and SIM. 



3 INTERFEROMETRY 

The opacity patch representing the non-uniform optical 
throughput of a telescope also has an impact on the per- 
formance of an interferometer. The orientation and phase 
distribution of the compressed beam is affected as derived 
in Sec. [2j and this will contribute to the error budget 
of visibility and phase determination, depending on im- 
plemen tation aspects of th e transport and combination 
optics (Buschcr et ahllioOSl ). Additionally, a specific con- 
tribution to baseline noise can be identified. 

An effective pupil centre position can be defined for 
each opacity patch, as average weighted by the transmis- 
sion distribution: 

_ / d£,dri ^ ■ tj^, V) _^ po Ap 

J d^dr, t{^,r!) To At -A,- 

where Ap and At are the patch and primary mirror areas. 

Thus, the pupil centre displacement is related to the 
position of the opacity patch, weighted by the relative 
opacity vs. unperturbed transmission, and by the relative 
size of the patch vs. pupil. 

The opacity patch centre position can range over 
the whole geometric extension of the pupil, with uniform 
distribution, so that the average pupil displacement is 
zero: 

(Cp) = J didn = , (6) 
and its RMS value is related to the RMS pupil size: 

(Cp) = (fe - {i.) f) = YoA^fi''^ ■ 

In simple geometry cases, the RMS size of the aper- 
ture can be easily expressed in terms of the linear size 
in the relevant direction, for a rectangular pupil, or 
of the diameter D, for a circular unobstructed pupil, i.e. 
respectively ^rms = Li^l^/Vl and ^rms = D/A. 

Non-uniform variations of the optical throughput, 
therefore, induce apparent changes in the position of each 
telescope of an interferometric array. 



3.0.1 PRIMA 

The baseline is of order of 100 m, and the stellar sep- 
arations considered for high precision astrometry in K 
band (2.2 ^m) are up to a few tens of arcsec. We con- 
sider a measurement performed by means of the Auxil- 
iary Telescopes (AT, diameter 1.8 m), and assume one 
opacity patch with 1% throughput degradation and size 
o = 5 cm, in an arbitrary position over the pupil. Then, 
the RMS baseline error is about 3.5 mm, and the cor- 
responding error for baseline B — 100 m and separation 
if = 10" is o-{ip) = 35 fias. 

3.0.2 SIM 

The baseline is _B = 6 m and the operating wavelength 
is in the visible, so that we select A = 600 nm; the indi- 
vidual apertures have diameter D = 0.3 m. Narrow angle 
astrometry is performed on separations up to ^ = 1°, 
so that, assuming an opacity patch with 1% through- 
put degradation and size a = 1 cm, we get a RMS base- 
line error of 0.7 mm and angular measurement error of 
o{ip) = 0".4. 

The larger astrometric error with respect to the 
PRIMA VLTI case is due to the combined contributions 
of shorter baseline and larger separation angle. 

3.1 Calibration vs. metrology 

The above values of error introduced by non-uniform 
transmission variation are not immediately related to the 
final measurement performance, because diagnostics and 
correction are not yet considered. Of course, on-sky cal- 
ibration procedures are foreseen, both for PRIMA and 
for SIM, in which the instrument is rapidly switched be- 
tween the science target and suitable reference objects. 
The errors discussed above are systematic and strictly 
correlated (although not equal) among different sources, 
and this is a key factor in the definition of a measure- 
ment sequence leading to their suppression. Depending 
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on availability of specific calibrators (e.g. binary systems 
or resolved objects considered sufficiently stable), simple 
calibration sequences can be devised. 

For example, if ip\ and (p2 are the angular quantities 
of interest referred to the science and reference sources, 
respectively, the measurement of the latter changes be- 
tween epochs ti , t2 only because of instrument variation 
by (7(992), i.e. 1/32(^2) = V'2(ii)+o"((y32), whereas the former 
changes by both astrophysical and instrument amounts, 
respectively (5(^1 and(7((^i): V5i(t2) = 'pi{t-C)+5ipi+a{ipi). 
From Equation |5] instrumental errors are correlated: 

a(^i) =^(^2) — , (9) 
<y32 

SO that the set of measurements allows correct determi- 
nation of the desired astrophysical quantity Stpi. 

In practice, the calibration procedure must be more 
complex, because the systematic errors affecting the mea- 
surement are likely to include additive components, as 
well as multiplicative components as from Equation [H] 

It should be noted that a metrology system installed 
to monitor the distance among fiducial points on the tele- 
scope structure is not sensitive to effective pupil displace- 
ments induced by non-uniform throughput variations, be- 
cause the change involves optical parameters and not the 
instrument geometry. Metrology can be extremely useful 
with respect to short term disturbances from the envi- 
ronment, and to long term geometry perturbations, but 
is liable to failure in the detection of the effects discussed 
herein. 



4 DISCUSSION 

The astrometric effects induced by non-uniform varia- 
tions of the telescope throughput depend of course on 
the characteristic size and amplitude of the defects on the 
mirror surfaces. In turn, they depend on details of man- 
ufacturing and operating environment which are not eas- 
ily identified nor measured. Useful practical indications 
may be provided by testing representative components 
e.g. during or after a significant period passed within an 
environment comparable to that of operation, as a cryo- 
genic vacuum cha mber, and subject to t he expected ra- 
diation doses, as in lBaccaro et al.l (l2004h . 

4.1 Design and implementation aspects 

For future astrometric instruments, using a comparably 
large field of view, the balancing of errors due to sym- 
metric aberrations could be used to advantage, allow- 
ing compensation over a set of measurements taken in 
suitable conditions (i.e. cases 1 and 4 above). An aber- 
ration control approach could be conveniently adopted 
at the early design stage, minimising the sensitivity to 
manufacture tolerancing and alignment errors, also with 
particular regard to the optimisation criteria of in-orbit 
re-alignment, to preserve the symmetry. From a strictly 
experimental standpoint, it is possible to devise concep- 
tual schemes for throughput measurement devices to be 
used on ground, for performance verification during the 
instrument integration, and even in orbit, to monitor the 



in-flight response. Their cost and benefit with respect 
to purely astronomical calibration is a matter of project 
trade-off. 



4.2 Impact on Gaia 

For Gaia, a significant fraction of the astrometric ef- 
fect associated with the image shape evolution induced 
by telescope throughput non-uniform variation could be 
taken care of within the framework of the self-calibration 
process, assuming that the instrument variation is suffi- 
ciently small over a time scale comparable with the cov- 
erage of the whole sky or a large part of it, e.g. order 
of six months. Also, part of the error may be averaged 
down among subsequent measurements, thus resulting in 
partial compensation. 

Partial compensation can be achieved also in case of 
sy mmetric aberrat i ons o ver the field of view, as discussed 
in lBusonero et all (|2006l ) on the subject of chromatic as- 
trometric errors. If opposite regions on the focal plane 
are affected by aberrations with opposite sign (as in cases 
1 and 4 above), the astrometric errors introduced on the 
two measurements cancel out. Also in case of partial sym- 
metry, or of limited balancing of the aberrations over the 
focal plane, some compensation can be achieved. 

Since the aberration distribution (and image profile) 
change over the focal plane, independently for each tele- 
scope, non-uniform throughput variations introduce as- 
trometric errors in the large angle measurement. The dis- 
tribution of photo-centre displacements in case of differ- 
ent WFE, as in figure O can provide an indication of the 
effect amplitude, but it is not fully representative because 
the patch positions are not correlated between the tele- 
scopes. The result is a field dependent error on the base 
angle, which is not strictly correlated to the measure- 
ments from the Base Angle Monitoring device (BAM), 
which probes only a small region of each optical system. 

An indication of throughput variation could be 
achieved by evaluation of the photometric response over 
a significant set of bright stars during each period of ob- 
servation, in order to monitor the instrument response 
at the sub-milli-mag level (corresponding to 0.01% over- 
all transmission variation). However, photometry is not 
sufficient per se to provide adequate information, since it 
cannot distinguish the cases of uniform and non-uniform 
variation over the aperture. 

A wavefront sensor (WFS) of the Hartman-Shack 
type, splitting the coUimated beam in several smaller re- 
gions, is potentially able to provide information on the 
variation of instrument transmission over its pupil, by 
photometric monitoring of the intensity from each region, 
as well as WFE derivative diagnostics by the elemen- 
tary image displacements on the focal plane. However, a 
WFS is sensitive also to the variations of the non-common 
part of the system, namely the geometry and response of 
the lenslet array. As for photometry based diagnostics, 
the detector might also suffer changes which should be 
identified and factored out from the estimate of optical 
throughput variation. 

Besides, since the astrometric effect is associated 
with a variation of the effective signal profile, suitable 
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image descriptors may be monitored, which usually have 
magnitude dependent sensitivity (e.g. image moments). 
The develop ment of diagnostic tools based on image pro- 
file variation ijGai an d Cancollicrc 2005) might prove use- 
ful to identification of critical instrument changes, in sup- 
port to the previous approaches. 

Future investigations might tackle the development 
of efficient monitoring and correction algorithms, based 
on the available information, as well as the investigation 
of cumulative and field-dependent astrometric effects as- 
sociated with the evolution of the telescope throughput. 
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Perryman M.A.C., de Boer K.S., Gilmore G., E[0g E., 
Lattanzi M.G., Lindegren L., Luri X., Mignard F., Pace 
O., de Zeeuw P.T., 2001, A&A, 369, 339 
Unwin S.C., Shao M., Tanner A.M. et al., 2008, PASP, 
863, 38 

Vucina T., Boccas M., Araya C., Ahhee C., 2006, Proc. 
SPIE, 6273, 62730W 



5 CONCLUSION 

Non-uniform variations of the optical throughput over 
the telescope aperture, as experienced in ground based 
instruments, may provide astrometric effects significant 
at the level of several microarcseconds. This may in turn 
represent a non negligible contribution to the error bud- 
get of the most challenging astrometric experiments cur- 
rently under implementation or proposed for the near 
future, as Gaia, PRIMA and SIM. 

The scale of the effects depends on manufacturing 
and operation parameters not easily measurable at the 
relevant level, i.e. order of 1% on mirror reflectivity over 
regions of few cm, or 0.02% on global throughput. Pro- 
cedures of periodic calibration on sky can be used to re- 
tain acceptable levels of the systematic errors introduced 
by instrument response variations, in particular for long 
term measurements as interferometric orbits. 

For Gaia, the time scale of instrument evolution must 
be compared with the period of effective parameter cal- 
ibration, to verify that the systematic errors are ade- 
quately removed from the data. Information from the 
WFS and from photometry might be exploited to moni- 
tor astrometric effects, as well as diagnostics of the image 
profile evolution. 
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